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ABSTRACT 



To elucidate the influence of shear flow on the generation of magnetic field (the 
dynamo effect), we study the kinematic limit where the magnetic field does not 
backreact on the velocity field. By non-perturbatively incorporating the effect 
of shear in a hehcally forced turbulence, we show that turbulence intensity and 
turbulent transport coefficients (turbulent viscosity, a and /3 effect) are enhanced 
by a weak shear while strongly suppressed for strong shear. In particular, /3 
is shown to be much more strongly suppressed than a effect. We discuss its 
important implications for dynamo efficiency, i.e. on the scaling of the dynamo 
number with differential rotation. 

Subject headings: MHD - stars: magnetic fields - turbulence 
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It is now widely accepted that astrophysical and geophysical magnetic fields are not 
the remains of a fossil field created during the formation of planets or stars (as they would 
have decayed on a time-scale much shorter than their current life time), but are self-excited 
by motions of conductive fiuid (for instance, molten iron within the outer liquid core for 
the Earth and conducting plasma for the Sun). The evolution of a magnetic field B in a 
conducting fiuid V is governed by the induction equation: 

9tB + V-VB = B- VV + r?V^B and V ■ B = , (1) 

where t] is the ohmic diffusivity. The first term on the right hand side (RHS) of Eq. ([1]) is 
the stretching of magnetic field lines by gradients of the velocity field. 

While laminar fiows that can generate magnetic fields (dynamo) have been known for 
a long time, the effect of turbulence on the generation of large-scale coherent magnetic 
field remains controversial. A main problem is that turbulence tends to create magnetic 
field at small scales (i.e. scale comparable to the original velocity field) while observations 
of astrophysical magnetic fields (for instance galaxies) reveal coherent magnetic field on a 
scale much larger than the fiuctuating velocity field. Theories, such as mean-field dynamo 



(IMoffatt 



1978 



Krause fc Radler 



19801), have investigated the necessary ingredients for 
large-scale field generation. In the framework of mean-field dynamo, the magnetic and 
velocity fields can be decomposed into mean and fiuctuating parts: B = (B) + b and 
V = (V) + V, where the (•) stands for an average on the realization of the small-scale fields. 
Substitution of this decomposition into Eq. ([T]) and averaging yield the following equation 
for the mean magnetic field: 

^^(B) + (V) ■ V(B) = (B) ■ V (V) + r]V^{B) + VxS . (2) 

The first term on the RHS of Eq. ([1]) represents the stretching of magnetic field lines by 
gradient of the mean fiow (V (V)) and is called the Q effect. It is an efficient mechanism 
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to cr eate toroidal field from a poloidal field in a system with differential rotation (IMoffatt 



19781 ). The term £^ = (v x b) is the electromotive force, which is often taken to be linear in 



the mean magnetic field ((B)). In the case of an isotropic turbulence, this can be simplified 



as: 



= a(5,)-/3(Vx (B)), + ... 



(3) 



The structure of the electromotive force permits the possibility of other mechanisms for 
the amplification of the large-scale magnetic field besides the Q effect. The one that has 
been discussed most is the a effect, the first term on the RHS of Eq. ([3]). This has been 
shown to generate magnetic field at large scale for a helical turbulence. Thus, it is a perfect 
candidate to explain magnetic fields in systems influenced by Coriolis force (which produces 
a net helicity) such as in stellar convection zones. This type of dynamo is thus classified as 
aQ if the Q effect (measured by the strength of the shear Q in our notations) is stronger 
than the a effect, or type if the a effect dominates over the fl effect. The second term in 
the RHS is Eq. ([3]) is the turbulent diffusivity which adds up to the molecular diffusivity rj. 
Consequently, if /3 is positive, it inhibits the growth of magnetic field. 

Recently, numerical simulations have shown dynam o action at large scale in non-helical 



turbulence in the presence of shear (lYousef et al. 



20081 ). This is an interesting result 



as the a effect is often thought to vanish in a turbulence without helicity. Various 
mechanisms h ave been invoked to explain this large-scale dynamo: stochas tic a effect 



(IProctor 



20071 ). shear amplification o f small-scale dynamo (jBlackmanlll998l ). magn etic 



effect driven by current helicity flux (iVishniac fc Cho 



20011) or negative diffusivity 



Urpin 



20031) which 



20021 ) . Another possibility is the shear current effect (iRogachevskii &: KleeorinI 
appears in a turbulent flow with a mean shear flow. In that case, the expression of the 
P coefficient can be rewritten Pijk = —0^eijk + Fjjfc(VUo) where (3'^ is the turbulent 
magnetic diffusion while the second term propor tional to shear VUq acts as a source 



of magnetic field (IRogachevskii fc Kleeorin 



20031 ). It is thus of prime importance to 
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investigate how the electro motive force (and consequently t 



on a large-scale shear flo w (iRogachevskii fc Kleeorin 



Brandenburg et al. 



2003 



le Qi and coefiicients) depend s 



2004 



Radler fc Stepanov 



2006 



20081 ) ■ In all these previous studies, strong shear is conductive to 



dynamo as it creates magnetic energy via the Q effe ct, acts as a source of m agnetic fleld 



(e.g. via the shear-current effect), causes instabihty (jTobias fc Hughes 



200J), etc. 



One interesting problem, which has not been investigated by most previous authors, 
is the effect of a stable shear flow on turbulent transport through the modiflcation of 
the properties of turbulence alone, without direct influence on (B) (i.e. no fi-effect, 
shear-current effect). A strong shear flow, without altering (B) direc tly, can reduc e 



19971) • This 



turbulent transport as turbulence becomes weak by shear stabilization (iBurrelll 
is basically because shear advects turbulent eddies differentially, elongating and distorting 
their shapes, thereby rapidly generating small scales which are ultimately disrupted by 
molecular dissipation on small scales (see Fig. [1]). As a result, turbulence level as well as 
turbulent trans port of various quantities can be signiflca ntly reduced compared to the case 



without shear (iKim 



2005 



2006 



Leprovost fc Kim 



20061 ). In particular, in the case when 



a stable shear flow is parallel to the magnetic fleld, a dramatic quenching of turbulent 
magnetic di ffusion (/3-effect) was clearly shown in a recent numerical simulation of 2D MHD 



turbulence ( iNewton fc Kim 



20081 ). In 3D MHD tur bulence, by considering a stable shear 



Leprovost &: Kim 



(120081 ) theoretically 



flow parallel to a uniform large-scale magnetic fleld, 
predicted that the a effect is quenched by shear as well as magnetic fleld. In particular, 
in the kinematic case (for weak magnetic fleld), the a effect was shown to be reduced as 
flow shear A increases with the scaling A^^^^. However, to understand fully the effect of 
shear on the dynamo process, it remains to compute its effect on the turbulence diffusion of 
magnetic fleld, i.e. the jS effect, by considering a non-uniform magnetic fleld. This is what 
we do in the remainder of this letter. 
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In the kinematic limit, the backreaction of the magnetic field on the velocity is 
neglected. From the physical point of view, this amounts to considering a very weak 
magnetic field and ignoring the Lorentz Force on the fiuid which is quadratic in the magnetic 
field. For an incompressible conducting fiuid, the resulting equations of motion are: 

9tV + V-VV = -Vj9 + z/AV + f, (4) 
dtB + y = B ■ VV + r]AB , 

V- V = V- B = . 

Here B is the Alfven speed, p is the total (hydrodynamical + magnetic) pressure and f is 
a small-scale forcing. To study the effect of shear flows and magnetic fields on small-scale 
turbulence, we prescribe a large scale fiow of the form Uq = —xAsy and a sheared 
large-scale magnetic field Bq = [Bq — Bx)ey. Bq has been chosen parallel to Uq so that 
there is no direct interaction between th e two fields, e.g. excluding the r2-effect in ou r study 



(in contrast with the case considered by (lYousef et al. 



2008 



Schekochihin et al. 



20081), etc.). 



To solve the equations for the fluctuating velocity fleld, u = V — Uq, and magnetic fleld, 
b = B — Bq, we use the quasi-linear approximation assuming that the interaction between 
fluctuating flelds is negligible compared to the interaction between large and small-scale 
fields. The equations for the fiuctuating fields can then be written as: 

dtn + u- VUo = -Vp + z/Au + f , (5) 
dth + u- VBo + Uo • Vb = b ■ VUq + Bq ■ Vu + r]Ah , 
V ■ u = V ■ b = . 

In the sequel, we shall assume an unit magne tic Prandtl number {u = rf) and introduce a 



time dependent Fourier transform (iKim 



20051) 
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Transforming the time variable from t to r = kx{t)/ky = kx{to)/ky + A{t — to), Eq. ([5]) can 
be written: 

drUi - Ux Si2 = -ikyOip - i{g^ + T'^)ui + fi , (6) 
drbi - TZuxSi2 = -hx6i2 + TZdrUi + ^7M^ - ^{g^ + r'^ybi , 
TUx + Uy + [3uz = rhx + hy + [3hz = . 

Here, TZ = B/A and 7 = B^ky/A are the ratio of the magnetic shear and constant magnetic 
field to the velocity shear, respectively; /? = k^/ky and g"^ = 1 + ^ = uk^/A and 
6*4 = (r, 1,/?). Note that since the first equation of does not involve the magnetic field. 



the s olutio r i to V j is the same as in the hydrodynamical case ([Kim 



from 



2OO5I ) ■ Using the velocity 



KimI (120051 ). the magnetic fluctuations can be obtained from the second equation of 



(ED as: 



TO 



A 



n 



+ ^7{T(r)-^(^)}-7^e(r-^) 



dt 



TO 



-a' 



,G{t,T) 



[n{l - i{Q{T) - Q{t)}) + i^{t - t)] 



Here, 



(7) 



G(t,r) 



-e{Q(r)-g(t)}, g(x)=/x + xV3 



^2 + ^2 g2^ f2 

I(t, x) dx , and J2 = 



T{r) - Tit) 

I{t, x){g^ + x^) dx 



where T{x) = axctaii{x / g) / g . by can be obtained using incompressibility: by = —rbx — f3bz- 

Our main interest is in the total stress and the electromotive force, which determine the 
growth/decay of the large-scale velocity field and the large-scale magnetic field, respectively. 
First, the stress is S = {uxUy) — (bxby). This total stress consisting of Reynolds stress {uxUy) 
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and Maxwell stress (bxby) gives a turbulent viscosity z/^ in Navier-Stokes equation for 
large-scale flows, which enhances the molecular viscosity to u + u'^. For the assumed shear 
flow Uq = —Ax, the turbulent viscosity is given hj S = u^A. Second, for the magnetic field 
considered here, the electromotive force reduces to: 



= {uybz - ujjy) = axyBo , (9) 

Sy = {u^bx - Uxbz) = ayyBo , 

= {uxby - Uybx) = -(3B . 
Note here that only three coefficients a^,, and (5 are non- vanishing in our configuration. 



In part icular, phenomena such a.s the O x J ( iRadler &: Stepanov 



effects (IRogachevskii &: Kleeorin 



20061 ) and shear current 



20031 ). which have been advocated to generate magnetic 



field for non- helical turbulence subject to rotation and shear as noted previously, are absent 



here (in contrast with the simulations of lYousef et al. 



20081 ). Note that a shear-current 



effect could be studied by using a similar analysis but assuming the large-scale magnetic 
field to depend on z rather than x, which will be addressed in a future contribution. 

To calculate the correlation functions involved in the transport coefficients, we consider 
an incompressible forcing which is spatially homogeneous and temporally short correlated 
with the correlation time r/. Specifically, in Fourier space, the correlation function of the 
forcing is taken as: 



(/^(ki,ti)/,(k2,t2)) = r/ (27r)='5(ki + ka) 6{h - h)(t>iA^2) 



(10) 



where the tilde denotes a Fourier-transform with respect to the spatial variable. As noted 
previously, the a effect can be linked to the helicity of the turbulent flow. Consequently, we 
consider a forcing with both a symmetric part (with energy spectrum E) and a helical part 
(with helicity spectrum H) given by: 



km{^) = E{k) (^5i^ - + ieimpkpH 



{k). 



(11) 
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In the following, the turbulence intensity, turbulent viscosity and a effect are expressed in 
terms of their values in the absence of shear or magnetic field, Cq, vi^, and /3o, which can 
be shown to be: 



eo 
/So 







Using equa tions for velocity in 



(2vr)2 

^/ 
(27r)2 

(27r)2 



V 

m 



(12) 



dk 



m 



KimI (120051 ) and Eq. ([71) and after a long algebra following 



Kinj (120051 ). we can find the turbulent intensity, stress and the electromotive force. Omitting 
the details, here we provide the results only for the hmiting case of a weak (^ = uk^/A ^ 1) 
and strong shear (^ = uky/A -C 1). 

First, in the case where the shear is weak compared to the diffusion rate (.^ ^ 1), we 
obtain: 



{hi) 



a 



xy 



2eo 
3 

eo 

eo 
3 

eo 
3 



1 + 
1 + 



9C 



35 
70 



(13) 



2 £—2 



+ 



2 , / S* 



35 

25267^2^-2 



715 



1 + 



4C: 

21 



a. 



yy 



P 



ao 
-l3o 



1 + 



33e 



-2" 



1 + 



70 
35 
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Note that the turbulent viscosity ut and the (3 effect are proportional only to the energy 
part of the forcing while the a effect is proportional only to the non-reflectionally symmetric 
part of the forcing. This is consistent with the expectation that the a effect is due to helical 
flow, which results from the helical forcing with helicity spectrum H. Eq. (fT3|) shows that 
(in the weak shear limit) all the turbulent coefficients increase with shear above their values 
without shear. The increase in / 3 with shear seems to be in agreement with numerical results 



shown in Fig. 1 of 



Mitra et al 



(120091 ) obtained in a slightly different configuration of Uq 
and Bq. Eq. f fT3|) also shows that a^y -C ayy i.e. that the electromotive force is primarily 
parallel to the large-scale magnetic field (i.e. in the y direction). Furthermore, without 
shear (^^^ = 0), we see that axy = showing that this component of the a effect exists only 
for non vanishing shear. This is due to the fact that shear induces an anisotropic turbulence 



[see e.g. 



Leprovost fc Kim 



20071 ) which in turn trigg ers off-diagonal compon ents in the a 



Kim &: DubruUd ( 120011 ) who found in 



tensor. Note that a different result was obtained by 
two dimensions th at the turbulent d i ffusiv ity decreases with shear. This difference comes 
form the fact that iKim fc DubruUd (120011 ) considered an anisotropic forcing, physically 



different form the isotropic forcing considered here. 

In the opposite limit of strong shear (^ = uk^/A ^ 1), turbulence intensity and 
transport coefficients are obtained as follows: 















ibl) 










Oixy 




^yy 





(14) 



These results show that in the limit of strong shear (compared to diffusion), all the turbulent 
quantities are reduced by shear with scalings given above. Note that the magnetic energy 
(6^) is more reduced than kinetic energy (m^). Furthermore, both the velocity and magnetic 
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field in the direction of tlie sliear are reduced more severely tlian in tlie perpendicular 
direction, manifesting the anisotropic turbulence induced by shear. It is because flow shear 
directly influences the component parallel to itself (i.e. the x component in Fig. 1) via 
elongation while only indirectly the other two components (i.e. the y and z components in 
Fig. 1) through enhanced dissipation. The electromotive force shows that the x-component 
of the a effect (a^j/) is now larger than the y one {oiyy). This is again because, as the shear 
increases, the anisotropy in the flow increases enhancing the off-diagonal component a^y 
strongly. Finally, the turbulent diffusivity (3 is reduced as ^"^^'^ more severely than the a 
effect {uyy oc i^^^)-, which has interesting implications for the dependence of the dynamo 
number (characterizing the efficiency of the dynamo) with differential rotation, as discussed 
in the introduction. 

To summarize, we found that the /5 effect is reduced as A~'^/'^, with a much stronger 
dependence on the shear than the a effect {uyy oc A~^/^). This result can have interesting 
implications for solar dynamo which is often envisioned to take place at the base of the 



convection zone where the shear is quite str ong (the so-ca 



for the weakness of the interface dynamo (jPikpati et al. 



ed_tachocline), e.g. to compensate 



20051 ). In particular, quenching 



by shear should be incorporated when assessing the efficiency of dynamo, e.g. the dynamo 
number given hy D = aQL^/{ri + /9)^, where Q is the differential rotation (corresponding 
to flow shear: Q = A in this paper) and L is a characteristic scale of the system. While 
i t is convention ally thought that the dynamo efficiency increases proportionally to shear 



( iKulsrud 



19991 ) for an af2 dynamo, our result suggests that the relation between the dynamo 
efficiency and the shearing rate is unlikely to be so simple. For instance, in the case of the 
aQ dynamo, the dynamo number D becomes: 

D = aALy{r] + f3f (xA\ (15) 

which increases with shear much faster than what has been conventionally thought. In the 



- 12 - 



case of an o? dynamo, we obtain a different scaling: 



(16) 



In the case of a^-dynamo, it is also interesting to examine how the growth rate of the 
magnetic field scales with shear: using standard formula for the maximum growth rate (see 



Moffatt 



1978 



for instance), we obtain the estimate o oc o? j oc ^ 



It is interestin g to note that our result s are y ery different from the recent works by 



Yousef et al. 



fl2008h and 



Schekochihin et al. 



(120081 ) where the large-scale magnetic field is 



amplified with a growth-rate scaling as J? or A. This is because, in these works, the 
dynamo instability is triggered by direct interaction between the large scale magnetic field 
(with both components parallel and perpendicular to the velocity field) and velocity field 
(i.e. (B) ■ V(U) 7^ 0). Note that in these works, the shear flow is assumed to be weak 
compared to the diffusion rate, corresponding to our weak shear limit (^ ^ 1). It would be 
interesting to study the opposite limit of a strong shear (^ <^ 1). 

Finally, we showed that turbulence and transport are enhanced for weak shear while 
quenched for strong shear. Therefore, there is a critical value of the shear for which the 



turbulence intensity and transport are maximum. As shown by lNewton fc Kind (120071 ). this 



can be due to resonance between the turbulence and shear flow when the characteristic 
frequency of turbulence matches the advection by shear flow (i.e. the Doppler shifted 
frequency vanishes). 
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Typical distance an eddy can transport a passive scalar field 



Fig. 1. — Sketch of the effect of shear on a turbulent eddy. 
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